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NUMERICAL INSTABILITY AND STABILIZATION OB THE KINEMATIC 
DIFFERENTIAL EQUATION USING QUATERNIONS 

W. Geiger 

Institute B for Mechanics, University of Stuttgart 
1 . Introduction /T 118* 

Because of the singularity of the Eulerian angles (frame barrier) 
other coordinates often must be used to describe the position of an, 
rigid body or '.gyroscope with a given rotating vector o> =Pqj(i) < respectively. 
One possibility is using rotating quaternions. These satisfy the 
linear differential equation (DGL) 
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The^jHare here the normalized coordinates of the rotating vector («y 
in the rigid system. The present text is concerned with this 
kinematic differential equation using quaternions ( KQD ) . 

The four coordinates of must fulfill the requirement 

i N - 1 = 2o + 2? j; Sl + $1 ~ 1 = ° I ( 2 ) • 

From KQD follows N = 0, thus N = constant is a first integral of 
KQD^ ^ . It is physically stabilized on the boundary. 

2. Numerical Instability of the Normalized Integral 

An unsymmetrical force-free gyroscope was employed as an 
example. Its motion is represented exactly by the Eulerian 
angles •. This concerns a periodic nutation ' in of the 


* Numbers in the margin indicate pagination in the foreign text. 

1 ) 

'The designation "integral" is also applied in the following for N. 
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figure axis. The Renge-Eutta procedure of the fourth order with 
constant step H was employed as integra.tion procedure. 



The normalized integral deviated in a nearly linear manner from 
the desired value. The course (designated hy KQD) is given as a 
function of time, where the first 30 periods of GEi were calculated 
with 30 steps for each period. Because of the physical stability 
of the normalized integral it is only a matter of a numerical in- 
stability here. The procedural error must always have an effect 
in one direction. As proof the quaternion from the disrupted 

Taylor series may be calculated after step (j+1) under the assump- 
tion (VI ) that the rotating vector may be assumed to be constant 
during an integration step with the result 

| II Vj+iII" = (1 *- p f ) \\rji\i 2 , v/ith, Pi = /i — (#IWjl)^V 

* I 8 k (3). 

t W , j = 0,1,2, ... ;• 


where for all possible steps p f >0. Corresponding values apply for 

procedures of the third and fifth order. The quaternion norm there- 
/ 

by decreases with each integration step and asymtotically approaches 
z ero . ~ ’ 

Remarks : 1151 

a) A rather good approximation of precondition VI was achieved', 
in the example . employed. The value of rotating vector could 
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even further he assumed as always constant (72). 

Then 

'^-‘const.i-Tf, ' ||7> = g (4) 

and this function describes the actual course of the normal- 
ized integral with a high degree of accuracy. 

b) If as an especially simple example a permanent rotation is 
calculated with ; Wi_= ‘a^o, !,,< ‘ equations (5) and (4) describe 
exactly the behavior of the normalized integral. 

c) The proof resulting from (3) requires in addition to 71 

only that the matrix nl of KQD is symmetrical and involu- 
2 ) 

toric ' , so that other differential equations with such 
coefficient matrices will exhibit identical characteristics. 

i ~ 

When the Euler i an angled is calculated from the quaternion 
according to the usual formula 

1 ws 0~ql~ qj-q * ( 5 ) , 

then the normalized error effects a very large error in this angle, 

for _ - 

i cos Hortp cos -i?* , ihus : $ » — 1 - a " 

— 1 - ■ tan . 

.... -- - ( 6 ) 

I — ^ 

when; is the exact angle. This causes a numerical singularity in 
the range of the frame barrier. 


Independent of these results, suitable steps must be taken to 
stabilize this integral numerically because of the numerical in- 
stability of the normalized integral. 


3 . Stabilisation according to Baumgarte 


Baumgarte [ 2 J has proposed a gradient stabilisation procedure 
for systems of usual differential equations of the first order. In 
the case of the requirement for stabilisation designated here "on 
average" the differential equation is altered with the aid of known 
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first integrals in such a manner that the stabilized differential 
equation leads to a differential equation for the error in integrals 
with the result that these are reduced asymtotically on the average 
when they have deviated from desired value zero through procedural 
error. Applied to kinematic differential equation using quaternions, 
the stabilized differential equation (KQDS) becomes 

q = i(n-4: WE)q (7), 

where k is a positive, otherwise arbitrarily selected stabilisa- 
tion factor. The normalized integral would have to behave according 
to 

LjV — 4fcA T (l + N)_ (g) 

Actually, however, the typical course shown in Pig. 1 with k = 0.5 
as an example occurs (designated by KQDS). During the first inte- 
gration step a jump in the norm occurs. Thereafter the integral 
stabiliyes on an erroneous value, the error is therefore not re- 
duced asymptotically. 


Examination of a Renge-Kunge step demonstrates that the series 

(N. ) of the normalized error fulfills the corresponding law of form- 
J 

ation 

l N j± i _= N) - (i +N ) )j) j , p L = Kj; Ny (9) . 


The expression for /^p can no longer be stated explicitly. It is a 

I 

polynominal of the 80th order in . The coefficients are in each 
case polynominals in k with parameters H and j-W . Here Y1 is applied 
again. 



Figure 2: Jump in Norm 
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Since the absolute member of does not disappear, the value 
of 1L already differs clearly from zero (Pig. 2), although = 0. 
Stabilisation would probably have no influence, as long as the 
solution meets the requirement. Because erroneous aids, subsequently 
used in the averaging, are, however, intentionally applied to the 
integration procedure, a correct solution via stabilisation is made 
incorrect. 

She numerical examination of the convergental behavior of the 
sequence (N^) suceeds, with the application of the additional assump- 
tion T2 although this entails a great deal of work. It converges 

V 

toward a value {j^|or -demonstrates divergent behavior according to the 
value of the stabilisation factor k^. She result is -presented in 
Pig. 2 (logarithmic . scaling, step size and rotating vector as in 
Pig. 1). Shere are only two optimal stabilisation factors, in p ar- /.S 120 
ticular the zero positions of In order to maintain- f^«>] < 30 -c j 

it is necessary to select the stabilisation -factor with a minimum 
of three decimal places exactly equal to one of the two optimal 
stabilisation factors. 

She criticism pertaining to the stabilisation procedure^ 
described in [2j may be summarized as follows: a) the accuracy 
depends to a great extent on the stabilisation factor; b) in some 
cases stabilisation factors have to be determined very exactly. 

She search for these definitely entails a great deal of work and 
these are above all dependent on the problem. Shey undergo changes 
with the differential equation supplied at the beginning, commencing 
conditions, integration procedure and stepsise chosen; e) proof 
assumptions are questionable from the standpoint of numerics. She 
asymtotic differential equation applied for the prbof of the error 
in first integrals (here: (8)) describes the behavior of the inte- 
grals only when the integration procedure functions with no error. 

Shis differential equation, however, is only actually fulfilled 

^^eojis a zero place of function |Q only dependent on IT. 

^Sh£ regulation designated here as "individual" stabilisation 
leads to. the same qualitative results.. 
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in the framework of the integration error. Since the commencing 
condition also remains within this order, it can scarcely describe 
the numerical behavior of first integrals. 

4. Stabilisation through Normalisation 

It has been shown that the difficulties resulting from the 
numerical instability of the quaternion norm are avdided completely, 
fthen the quaternion is divided by its norm after each integration 
step. A proof of this procedure will be published separately. 
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